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ABSTRACT 



Aims. We present an analysis of the star formation properties of field galaxies within the local volume out to a recession 
velocity limit of 3000 km s _1 . 

Methods. A parent sample of 863 star-forming galaxies is used to calculate a £?-band luminosity function. This is then 
populated with star formation information from a subsample of 327 galaxies, for which we have Ha imaging, firstly by 
calibrating a relationship between galaxy B-band luminosity and star formation rate, and secondly by a Monte Carlo 
simulation of a representative sample of galaxies, in which star formation information is randomly sampled from the 
. observed subset. 

Results. The total star formation rate density of the local Universe is found to be between 0.016 and 
0.023 Mq yr _1 Mpc -3 with the uncertainties dominated by the internal extinction correction used in converting 
measured Ha fluxes to star formation rates. If our internally derived B-band luminosity function is replaced by one 
from the Sloan Digital Sky Survey blue sequence, the star formation rate densities are ~ 60% of the above values. 
We also calculate the contribution to the total star formation rate density from galaxies of different luminosities and 
Hubble T-types. The largest contribution comes from bright galaxies with Mb ^-20 mag, and the total contribution 
from galaxies fainter than Mb =-15.5 mag is less than 10%. Almost 60% of the star formation rate density comes from 
galaxies of types Sb, Sbc or Sc; 9% from galaxies earlier than Sb and 33% from galaxies later than Sc. Finally, 75 - 80% 
of the total star formation in the local Universe is shown to be occurring in disk regions, defined as being >1 kpc from 
f^*) ' the centres of galaxies. 

Conclusions. The star formation rate density estimates found here are consistent with values from the recent literature 
using a range of different star formation indicators. Even though they are numerous, dwarf galaxies contribute little to 
the star formation in the local Universe, and the bulk of the star formation takes place in L* spirals. 
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1. Introduction sample all star-forming spiral and dwarf galaxy types (SOa 

- Im), and the galaxies have recession velocities less than 
In spite of many studies of the star formation (SF) process 3000 km s -i A11 ga i axies were observed with the 1.0 me- 
in all types of galaxy, and of the variation of the total SF tre Jacobus Kapteyn Telescope (JKT), part of the Isaac 
rate (SFR) as a function of cosmological look-back time, Newton Group of Telescopes (ING) situated on La Palma 
many questions remain. In particular, a full census of the in the Canary Islands. The s election and the obs ervation of 
SF activity across all types of field galaxies, including the thc sample are discusse d in lJames et all d200l . hereafter 
faintest star- forming dwarfs, is yet to be completed. The p a p e r I 
Ha Galaxy Survey (HaGS) is a study of the SF proper- 
ties of galaxies in the local Universe, using fluxes in the Since HaGS includes a large number of dwarf galax- 
Ha line to determine the total rates and distributions of ies, which are only detectable at very small distances, the 
SF within the selected galaxies, as a contribution to just sample as it stands is not volume-limited. However, the 
such a census. The observations cover 334 galaxies, which galaxy selection was performed according to well-defined 

criteria in terms of galaxy apparent magnitudes, diameters 

Send offprint requests to: P. A. James and recession velocities, and so it is possible to calculate in- 

* Based on observations made with the Jacobus Kapteyn completeness corrections for galaxies of all types, from the 

Telescope operated on the island of La Palma by the Isaac fraction of the total surveyed volume in which they could 

Newton Group in the Spanish Observatorio del Roque de los lie and still satisfy the selection criteria. By applying such 

Muchachos of the Instituto de Astrofisica de Canarias corrections, volume density statistics, such as luminosity 
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functions, can be calculated. This is the primary aim of the 
present paper. This analysis will then enable the calculation 
of the average SFR density (in units of M@ yr _1 Mpc~ 3 ) 
due to all galaxies, and the same quantity subdivided into 
the contributions from different categories of galaxy. 

Thus the key scientific questions we will address are as 
follows: 

How does the total SFR density derived from the HaGS 
sample compare with other local (but generally larger- 
scale) estimates, and with estimates at intermediate and 
high redshifts (recently reviewed by Hop kins fc Beacoml 
12006ft ? 

What fraction of the SFR density is contributed by 
galaxies as a function of their morphological type, and as 
a function of their optical luminosity? 
What fraction of the SFR density occurs within the central 
bulge regions of spiral galaxies, and what fraction in disks? 

This approach is complementary to that adopted by 
other surveys of SF in the low redshift Universe (e.g. 
iGallego et al.|[l995tlBri nchmann ct al. 2004|) in that it sam- 
ples low-luminosity spiral and irregular galaxies more fully 
than any previous study. The SF indicator used (narrow- 
band CCD imaging in the Ha line) also allows good esti- 
mates of total SFRs without the need to correct for missed 
flux which occurs for slit- or fibre-based spectroscopy, and 
gives two-dimensional spatial information which is not 
available from such spectroscopic approaches. 

The structure of the current paper is as follows. After 
summarising the survey data in Section 2, we derive, in 
Section 3, the local (B-band) luminosity function from 
the parent sample of galaxies used for the survey. This 
will be compared to luminosity functions from the litera- 
ture derived from a comparable field sample of galaxies. In 
Section 4, the local SFR density will be calculated, first by 
a direct conversion of £>-band luminosity to SFR, and sec- 
ondly using a Monte Carlo method to extrapolate the prop- 
erties of the observed sample of galaxies to the whole par- 
ent population from which they were drawn. In Section 5, 
the contribution to the total SFR density from galaxies 
of different absolute magnitudes and types is investigated. 
The split of SF between the central 1 kpc in radius and 
disk regions is explored in Section 6. Section 7 contains 
a comparison of the total SFR density estimates derived 
here with literature values, and all results are summarised 
in Section 8. Finally, the appendix looks at how the mor- 
phological types used in the present analysis map onto the 
widely-discussed red and blue galaxy sequences. 

2. The galaxy sample selection 

The primary aim of HaGS was to determine the total SFR 
of all galaxies in the local Universe, and to analyse the con- 
tributions to this total from galaxies of different types and 
luminosities. In particular, it was considered very impor- 
tant to give detailed consideration to low-luminosity dwarf 
galaxies, which are numerically the dominant galaxy pop- 
ulation, but are often omitted from studies of this type. 
For example, the fainter galaxies are inevitably numeri- 
cally under-represented in any magnitude-limited sample. 
In this Section we describe the sample selection used to 
overcome these problems, which enables us to infer the to- 
tal SFR of all galaxies from our observations of 327 galaxies 
(the remaining seven galaxies were serendipitously observed 



objects and are thus omitted from the analysis presented 
here). 

The HaGS sample w as selected f rom the Uppsala 
Galaxy Catalogue (UGC; iNilsonl I1973T) . We assume the 
UGC is complete within its selection criteria, that is, all 
galaxies down to a limiting diameter of 1.0' at declinations 
greater than —2.5°. Selection bi ases in the UGC are bey ond 
the scope of this work (but see lde Jong fc Lacevl l2000 who 
find incompleteness of both high and low surface brightness 
galax ies in the similar ESO-Uppsala catalogue, iLaubertsl 
11982ft . 

The HaGS selection criteria are described in detail in 
Paper I, but can be summarised as follows: type SOa (T = 0) 
or later; diameters between 1.7 and 6.0 arcmin and recession 
velocities less than 3000 kms" 1 . A £?-band magnitude limit 
of 15.5 was adopted, corresponding to the limiting magni- 
tude beyond which the UGC becomes significantly incom- 
plete. A Galactic latitude limit of \b\ >20°was also imposed 
to minimise the effects of Galactic obscuration. In addition, 
to prevent the skewing of the results by the high density 
and special morphological characteristics of galaxies in the 
Virgo cluster, galaxies in regions centred on M 87 were ex- 
cluded. We investigate the effects of varying the size of the 
excluded region in section 3, but our adopted exclusion for 
the initial sample selection lead to a sample (which we refer 
to as the parent sample) of 863 galaxies. The observing list 
was trimmed by removing the most highly-inclined spiral 
galaxies (a/b >4.0) since these suffer from the worst in- 
ternal extinction effects and yield less information on the 
spatial distribution of SF than do face-on galaxies. The as- 
sumption here is that galaxy orientations are sufficiently 
random that no overall bias is introduced by this cut, i.e., 
the remaining 743 face-on galaxies are statistically repre- 
sentative of the full parent sample. 

Recession velocities are available on the NASA/IPAC 
Extragalactic Database (NED) for all UGC galaxies that 
satisfy the diameter, apparent magnitude and type require- 
ments for our parent sample. There is one frustrating excep- 
tion, UGC 4413, which is just large and bright enough for 
conclusion, but has no measured recession velocity on NED. 
In any case, UGC 4413 is highly inclined and and would be 
excluded from the pool of galaxies considered sufficiently 
face-on for potential observation. Overall, it can be consid- 
ered that there is no incompleteness in the current sample 
due to missing recession velocities. 

Within the time allocation, we observed a sub-sample 
of 327 galaxies from the parent sample, which sampled the 
full range of galaxy parameters in the parent sample. This 
gave us detailed SF parameters for each galaxy within this 
sub-sample. The analysis presented in this paper uses these 
parameters to supply estimates of the likely SF properties, 
and the scatter in these parameters, of galaxies as a func- 
tion of type and absolute magnitude. These distributions 
are then used to populate a simulated volume of galaxies 
with SF properties, where the relative numbers of galaxies 
of different luminosities are drawn from a -B-band luminos- 
ity function derived from our parent sample. Whilst this 
process is likely to be significantly in error for individual 
galaxies, it should give statistically reliable results when 
averaged over many tens or hundreds of galaxies. 

The apparent blue magnitudes quoted by the NED for 
each galaxy were used. The NED magnitudes are gener- 
ally taken fro m the Third Reference Cata logue of Bright 
Galaxies (RC3: lde Vaucouleurs et al"1ll991ft . To convert the 
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Fig. 1. The distribution of absolute blue magnitudes within 
the parent sample of 863 galaxies, as a function of distance 
in Mpc. The solid line corresponds to the apparent magni- 
tude limit of B =15.5 mag used in selecting the galaxies. 
Black/filled symbols denote the observed subsample, and 
red/open symbols the remaining unobserved galaxies. 

apparent magnitudes to absolute magnitudes, distances 
were calculated for each galaxy using a V i rgo-in fall cor- 
rection method based on that of ISchechterl (|1980f) and an 
asymptotic Hubble constant of 70 km s _1 Mpc -1 . Such 
distances are very uncertain for the closest galaxies, and, 
where available, distances from the literature based on ei- 
ther Cepheid variables or the apparent magnitude of the 
tip of the red giant branch were used in preference. Fig. [T] 
shows the distribution of absolute blue magnitudes (Mb) 
within the parent sample as a function of distance in Mpc, 
with the solid line showing the limiting magnitude adopted 
of B = 15.5 mag. The B magnitudes plotted here, and used 
in the following Section for the calculation of visibility vol- 
umes, are not corrected for Galactic extinction. However, 
the B magnitudes used subsequently, for the calculation of 
the galaxy luminosity function and the estimation of global 
SFRs, have corrections applied based on the methods of 
Schlcgcl, Finkbeiner & Davis (1998), which are also used 
to correct all Ha measurements (see Section 4). 

3. The i?-band luminosity function 

The first stage of the analysis is to use the parent galaxy 
sample to produce a i?-band luminosity function for galax- 
ies in the volume sampled by this survey. We will then use 
the representative subsample observed in Ha emission to 
populate this £?-band function with information on SFRs 
of individual galaxies, from which we can calculate the to- 
tal SFR in the local Universe. This two-stage process was 
preferred to the direct calculation of the Ha luminosity 
function from the 327 galaxies with Ha data because the 
latter sample is too small to constrain fully the form of the 
luminosity function. 

The luminosity function can be most directly calculated 
from a volume-limited catalogue containing a fair sample of 
all relevant galaxies over the range of intrinsic luminosities 
of interest. This is not practicable, however, for the study 
of field galaxies over a wide range of intrinsic luminosity, as 
is the case with the present study; the volume over which 



the faintest galaxies can be seen is so small that, far from 
containing a fair sample of the brightest galaxies, it may 
contain none of them at all. In this case, it is necessary to 
use a statistical approach, based on the assumption that the 
local Universe does not differ significantly in its galaxy con- 
tent and volume density from the Universe averaged over 
successively larger volumes. At some level, this assumption 
must be in error, since we live within the density enhance- 
ment of the Local Group, and close to the Virgo cluster and 
hence to the centre of the Local Supercluster (although the 
Virgo cluster galaxies themselves have been removed from 
the sample). This is likely to bias the calculated volume 
density to high values, but this is unavoidable and should 
be borne in mind when interpreting the results. 

We use a stati stical method ba sed on t he class i cal pr oce- 
dures outlined bv lSchmidtl (jl968l ) and bv lFeltenl (|1976h . to 
derive a B-band luminosity function from the parent sam- 
ple in our present study. Mathematically, this is done by 
numerically weighting each of the parent sample galaxies 
by (1/Vm), where V m is the size in Mpc 3 of the volume in 
which the galaxy could lie and still be included in our sam- 
ple, satisfying the apparent magnitude, diameter and reces- 
sion velocity criteria. Then the luminosity function <j>{Ms) 
is given by 

4>(Mb) is the space density of galaxies with blue magni- 
tude Mb, is the visibility volume within which the ith 
galaxy could lie and still be included in the sample, and the 
quantity f2/47r is the fraction of sky covered by the sample 
in question. In this case, £1 = 2irx 0.695, once the Virgo 
Cluster, Galactic plane and southern hemisphere regions 
have been removed. 

Each selection criterion will define a maximum volume 
in which a galaxy of given properties can be observed. 
The visibility volume for each galaxy is thus the volume 
which satisfies all t hese criteria (jDisnev fc Phillippsl . 119831 : 
iDavies et allll994h . There are three selection limits appli- 
cable to the parent sample (apparent magnitude, diameter 
and recession velocity), each of which leads to constraints 
on the distance range within which a galaxy must lie to 
be included in the sample. We adopt a limiting apparent 
magnitude, mu m , of 15.5 mag, and have excluded the small 
number of UGC galaxies fainter than this in defining the 
parent sample. A galaxy with absolute magnitude M must 
therefore lie within a distance in Mpc given by 

dmax.mag (Mp C ) = lO - 2 ^"^" 25 ) . (2) 

This equation applies to magnitudes uncorrected for 
Galactic and internal extinction. There is no bright appar- 
ent magnitude limit so d m in.mag = Mpc for all galax- 
ies. The lack of bright galaxies at distances up to 16 Mpc, 
clearly apparent in Fig.[TJ is due to the upper diameter limit 
on galaxy selection, explained in the following paragraph. 

The HaGS sample was selected to contain galaxies with 
apparent angular diameters between 1.7' and 6.0'. A galaxy 
with an apparent angular diameter D will be excluded from 
the sample if its distance, d, is low enough for the apparent 
diameter to be greater than 6.0', or if it is so far away that 
the apparent diameter falls below 1.7'. Thus the distance 
limits within which it can be observed are 
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Galaxies which are intrinsically small and faint will only 
be detected at the smallest distances due to the apparent 
magnitude and the minimum apparent diameter limits. The 
galaxies with the largest intrinsic sizes can only be observed 
at larger distances, due to the maximum angular diameter 
limit. 

The visibility volume V m is then defined by 



V m (Mpc 3 ) 



-in 

T 



(rl 3 

v max.lim 



d 3 ) 

min.lim/ ) 



(5) 



where d max ,iim is the smallest of d m ax,ma g , d maXjD and 
the distance corresponding to 3000 km s _1 in the direction 
of the galaxy. d m i n ,ii m is equal to d m i n ,D- 

The B-band luminosity function is then calculated by 
performing the sum in Equation [T] over all galaxies in the 
parent sample. The resulting function is shown in Fig. [5] 
The error bars contain t wo contr i bution s: 1 a Poisson errors 
taken from the tables of lGehreli (|1986f ). which are added in 
q uadrature to distance err ors, calculated using the method 
of Ide Jong &: Lacevl (|2000f) . Here we adopt distance errors 
of ±50% for all galaxies within 5 Mpc; ±30% for all galaxies 
between 5 and 10 Mpc; ±20% for all galaxies between 10 
and 15.6 Mpc; and ±15% for all galaxies beyond 15.6 Mpc. 

Luminosity functions of field galaxies are often found to 
be described well by a Schechter function (|Schechterl . |1976[ ) 
(see, e.g., Efstathiou, Ellis & Peterson 1988): 



cj>{L)dL = (f>*(L/L*) a ex V {-L/L*)d{L/L*), (6) 

This function decreases as a power law with increasing 
galaxy luminosity at faint magnitudes, but cuts off sharply 
for galaxies brighter than some characteristic magnitude, 
M*, corresponding to L*. a gives the slope of the luminos- 
ity function at the faint end and <f>* sets the overall normal- 
isation of galaxy density. 

The luminosity function derived from our parent sam- 
ple is shown as the points with error bars in Fig. [21(a) . A 
Schechter fit was made to this luminosity function, and is 
shown in Fig. HJa) as the dot-dashed line. The values of 
M* and <fi* in this fit are poorly constrained, due to the 
small number of galaxies brighter than the turnover at M* 
and the degeneracy between these parameters. Formally, 
the best-fit values are M* = -20.49 ±0.19, </>* = 0.00280 ± 
0.00062, and the faint-end slope a = -1.44±0.06. It should 
be noted, however, that a single Schechter function is not 
a good fit, due to the plateau in the observed function be- 
tween Mb = —18.0 and -19.5 mag. 

The solid line in Fig. [2ja) is an estimate of the B- 
band luminosity function for blue-sequence galaxies. This 
is derived from the Sloan Digital Sky Survey (SDSS), us- 
ing the New York University Value- Added Galaxy Catalog 
(NYU-VAGC; IBlanton et~aLll2005bl) . in p articular, the low- 
z galaxy sample (iBlanton et all l2005af) . This sample in- 
cludes galaxies with r < 17.8 mag in the redshift range 
1000 < cz < 15 000kms _1 (14 - 212 Mpc assuming 
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Fig. 2. (a) The luminosity function for the entire par- 
ent sample, types SOa - Im, plotted as points with error 
bars; the solid line shows the luminosity function from the 
SDSS blue sequence, (b) Luminosity functions including all 
galaxy types: SDSS (solid line); this study with all Virgo 
galaxies included (dashed line); and this study with galax- 
ies excluded out to 10°from the centre of the Virgo cluster 
(dot-dashed line). 

Hq = 70 km s _1 Mpc -1 ). Galaxies have been eyeballed 
to check photometry and redshifts have been corrected to 
distances using a local flow model for the nearest galaxies. 
K-corrections are applied to bring the photometry to z = 0. 
The NYU-VAGC absolute magnitudes were converted to B 
magnitudes using 



B = . g + 0.14 + 0.37(5 -r) 



(7) 



which was obtained by m atching Mi ll enniu m Galaxy 
Catalogue (MGC) galaxies (|Liske et all f2003h to SDSS 
galaxies and fitting a straight line to B — g versus g — r 
(using z < 0.1 galaxies). The SDSS luminosity function 
was then determined by summing, in absolute magnitude 
bins, 

Y 1 ^ (8) 

where V max is obtained from the NYU-VAGC and C s b is 
the surface-brightness completeness for each galaxy esti- 
mated using table 1 of IBlanton et all (|2005afh The lumi- 
nosity function is then effectively complete to a half-light 
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surface brightness of ~ 24 mag arcsec -2 in r. Counts and 
magnitudes were all converted to Hq = 70 km s _1 Mpc -1 . 

To determine the SDSS luminosity function for late type 
galaxies three cuts were considere d: (i) a colour cut se lect- 
ing blue-sequence galaxies as per iBaldrv et al.l (|200l : (ii) 
Sersic index less than 2.5 from the NYU-VAGC; (iii) con- 
centration index, i?90/i?50, less than 2.5 from the SDSS 
pipeline. These three definitions of the blue sequence yield 
largely consistent luminosity functions, and we adopted the 
first, hereafter referred to as the SDSS blue-sequence lu- 
minosity function, to compare with that from the HaGS 
parent sample. 

Overall there is a fair agreement in the shape and nor- 
malisation of the SDSS blue-sequence and HaGS luminosity 
functions, but there are two areas of significant disagree- 
ment. The first is that the HaGS luminosity function is 
somewhat higher around Mb = —19.0 to —20.5 mag, which 
is significant as these galaxies dominate the stellar mass and 
thus, plausibly, the SFR density of the overall galaxy pop- 
ulation. Secondly, the HaGS luminosity function shows a 
somewhat steeper faint-end rise, i.e., a larger dwarf-to-giant 
ratio than the SDSS function. The possible causes of these 
differences will now be investigated. 

One possibility is that the higher amplitude of our lumi- 
nosity function is due to the presence of the Virgo cluster, 
which may still have a significant influence even though the 
core region is excluded. This possibility was explored by 
producing luminosity functions for local galaxies selected 
from the UGC, with differently-sized regions of the Virgo 
cluster removed. For ease of comparison with literature 
studies, this comparison was done for all Hubble types, i.e. 
with types earlier than SOa reinstated, but otherwise the 
selection criteria are as described above. Two luminosity 
functions derived in this way are shown in Fig. (Hb), where 
the dashed line includes all Virgo cluster galaxies, and the 
dot-dashed line shows the effect of excluding all galaxies 
within a 10° radius centred on M87. (5°, 15° and 20° radius 
regions were also investigated, with similar results). For 
comparison, the full SDSS luminosity function, including 
all galaxy types, is shown as the solid line. The main effect 
of Virgo exclusion on the UGC luminosity function is to de- 
press it between Mb values of -18 and -20, thus causing the 
plateau mentioned above. The bright end of the function is 
then in fair agreement with the SDSS-derived luminosity 
function. At Mb = -20, the UGC-derived function has an 
a value </>(-20) of 0.0031 Mpc~ 3 mag~ 1 with Virgo included, 
falling to 0.0029, 0.0027, 0.0027 and 0.0025 Mpc^mag- 1 
for excluded regions of radius 5°, 10°, 15° and 20° re- 
spectively. For the full SDSS-derived luminosity function, 
</>(-20) is 0.0024 Mpc _3 mag _1 . This improved agreement 
shows that at least part of the bright-end discrepancy in 
Fig. [DJa) is due to the SDSS blue sequence definition ex- 
cluding a larger fraction of galaxies than our SOa - Im type 
constraint. This is supported by the analysis presented in 
the appendix, which shows that the blue sequence most 
nearly corresponds to types Sb or later. 

The literature sa mple most comparabl e to our par- 
ent sample is that of ISantiago et alj {l996), who studied 
the luminosity and diameter functions of galaxies within 
8000 km s _1 selected from both the UG C and its south - 
ern counterpart, the ESO-Uppsala Survey (|Laubertsl [l982). 
They do not present Schechter function fit parameters for 
their luminosity functions, citing significant biases in the 
sample selections which are hard to quantify, but the </>{— 



20) values for their luminosity functions from both cata- 
logues are somewhat lower than those we derive, at 0.0019- 
0.0021 Mpc _3 mag _1 , after correction from their assumed 
Hubble constant of 100 km s _1 Mpc -1 . 

It is clear from this analysis that cosmic variance may 
have a significant effect on the HaGS luminosity function 
over the volume we survey, and this will be more important 
at the faint end, due to the even smaller region sampled 
by the faintest galaxies. In order to quantify this faint-end 
difference, estimates were made of the amplitude of the B- 
band LF at Mb = -16.0, for comparison with literature 
studies. With all galaxy types and the Virgo cluster in- 
cluded, our UGC-based sample gave a value for </>(-16) of 
0.030 Mpc~ 3 mag _1 falling to 0.025 Mpc" 3 mag" 1 with the 
exclusion of the Virgo cluster, and to 0.020 Mpc~ 3 mag _1 
with the additional exclusion of ga l axies earlier than type 
SOa. The SDSS and ISantiago et al.l (|1996f) luminosity func- 
tions are significantly lower in this region, with 0(— 16) val- 
ues between 0.006 and 0.011 Mpc _3 mag _1 . Thus cosmic 
variance seems likely to have a substantial effect at these 
faint magnitudes. 

However, it is also true that the larger-scale surveys have 
substantial incompleteness below Mb ~ — 16 mag, and the 
corresponding uncertainties make it impossible to rule out 
the higher faint end counts that we find in the current sam- 
ple. Indeed, some other studies of the local galaxy luminos- 
ity function support a faint end slope as steep as, or steeper 
than, t hat found here for the HaGS parent sample. For ex- 
ample, [Marinonr^rrn] (|1999f ). studying an all-sky sample of 
6400 optically selected galaxies within 5500 km s _1 , found 
a very steeply-rising faint end slope (a < -2 for the faintest 
galaxies, vs ~-l.l over most of the luminosity range stud- 
ied) with the steepening at Mb ~ — 15 mag being ascribed 
to Magellanic irregular galaxies. 

In light of these uncertainties in the form of the local 
field galaxy luminosity function, for the remainder of this 
paper we will adopt the HaGS and SDSS blue sequence 
luminosity functions as limiting cases, and use both when 
characterising the SF properties of the local Universe. 

4. The local SFR density 

4.1. Methods for estimating the SFR density from the 
B-band luminosity function 

The B-band luminosity function calculated in the previ- 
ous Section from our parent sample gives the total number 
of galaxies of different luminosities per unit volume of the 
Universe. The next stage of this analysis is to populate the 
individual galaxies making up this luminosity function with 
SF information, thus enabling calculation of the total SFR 
density, and the breakdown of this total by galaxy type. 
The SF information is provided by the observed subsample 
of HaGS galaxies. 

Two methods are presented here. The first uses galaxy 
.B-band luminosity as a direct proxy for SFR, with calibra- 
tion factors linking the two quantities being derived from 
the observed subsample, and applied to the full parent sam- 
ple. The scatter in this conversion process is significant, 
however, and there are possible concerns that taking mean 
properties of galaxies may underestimate the effects of out- 
liers in the distributions of SF properties. These problems 
are addressed with the second method, which uses a Monte 
Carlo method to build multiple models of the Universe, 
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with galaxy types and luminosities drawn from the parent 
sample, and SF information randomly sampled from galax- 
ies with similar properties in the observed subsample. 

The Ha fl uxes are converte d to SFR values using the 
calibr ation of iKennicuttl (Q~998), which assumes a lSalpete"rl 
(|l955f) stellar initial mass function. 

4.2. The SFR:Lg correlation 

The first of these two methods relies on the fact that 
B-band luminosity (Lb) in star-forming galaxies largely 
comes from the young stellar population, and hence we can 
expect galaxy Lb values to correlate strongly with SFRs. 
This leads to a simple way to calculate the total SFR den- 
sity: we calibrate the SFR:Ls relation from the observed 
subsample, then use it to replace the Lb of each galaxy in 
the parent sample with a SFR. Finally, the same 1/V max 
analysis as was used to calculate the -B-band luminosity 
function will now result in the total SFR density. 

Figure [3] shows the log of galaxy SFRs per unit -B-band 
luminosity (in solar units), plotted against the log of the 
B-band luminosities, with different symbols showing early 
Hubble types (T-type 0-2; circled points), intermediate 
types (T-type 3-7; stars) and late types (T-type 8 - 10; 
filled circles) . All measured fluxes are corrected for Galactic 
extinction, and the SFRs have been corrected for extinction 
internal to the galaxy concern ed using the type-d ependent 
corrections derived in Paper II (I James et al.ll2005ft . Overall, 
Fig. [3] shows no strong trend, confirming that .B-band lu- 
minosity can be used as a proxy for SFR, but it is clear 
that, for example, the bright, early type galaxies have a 
lower mean ratio than is typical of other galaxies. Thus, 
calibration factors were calculated for each of these three 
type bins (T =0 - 2, 3 - 7 and 8 - 10) and three Lb bins 
(Faint, log(Ts)<8.5; Intermediate, 8.5<log(Ts)<9.5; and 
Bright, log(Ts)>9.5). In each case, the values were calcu- 
lated by dividing the sum of the galaxy SFRs by the sum of 
their B-band luminosities in solar units, and finally taking 
the log. The resulting values are given in Table [T] 

It is now simple to convert the Lb values for each 
of the parent sample galaxies to SFR, by multiplying by 
the calibration factor for the appropriate type and lu- 
minosity bin. Applying a 1/V m weighting and summing 
over the whole sample, we find the SFR density to be 
0.0237±0.0030 M yr" 1 Mpc~ 3 . Here the quoted error is 
the random error resulti ng from uncertaint ies in the indi- 
vidual galaxy Ha fluxes (j James et all 120041) and distances, 
and Poisson errors on the numbers of galaxies used in cal- 
culating mean binned properties. Much larger systematic 
uncertainties result from the extinction corrections applied 
in the conversion of Ha fluxes to SF rates, and from the 
luminosity function used. These factors will be addressed 
in the remainder of this section. 

We will next find the SFR density that results from 
using the SDSS blue sequence luminosity function. Again 
the method employed was to replace Lb with SFR values 
from our observed galaxy sample. However, in this case 
we only have the B-band luminosity function, and no in- 
formation on the fraction of galaxies of each Hubble type 
within each Lb bin. Therefore, we calculate the (SFR/T^) 
calibration factors simply as a function of Lb within our 
observed subsample. Since the sample is not subdivided by 
types, we can bin more finely in Lb, and so the mean cal- 
ibration factors are calculated within nine bins of width 




Log (B luminosity) 

Fig. 3. The log of the SFR to B— band luminosity ratio, 
plotted against the log of the B— band luminosity, for the 
observed galaxies. Circled points denote early types, T =0 
- 2; stars are intermediate types, T =3 - 7; and filled circles 
are late types, T =8 - 10. 



0.4 in log L B . The SFR density from the SDSS B-band 
luminosity function using this 'direct' calculation method 
0.0138±0.0015 M Q yr" 1 Mpc" 3 . 

This is significantly lower than the value from the 
internally-derived luminosity function, due to the overall 
lower level of the SDSS luminosity function. 



4.3. A Monte Carlo method for deriving the SFR from the 
B-band luminosity function 

There is a concern that the above method implicitly as- 
sumes a perfect correlation between B-band luminosity and 
current SFR; this is clearly not a valid assumption, and even 
within our type-luminosity bins, the scatter in (SFR/Lb) is 
large, at 0.3 dex. We have investigated the possible effects 
of this scatter using a Monte Carlo method, producing mul- 
tiple realisations of the local Universe. These models were 
populated using the 1/V mi T-type, and Mb information for 
each of the 863 galaxies in the parent sample. The 1/Vm 
value gives the number of times each galaxy is reproduced 
in the model universe. For each of these reproduced galax- 
ies, a SFR is selected by randomly sampling a galaxy from 
the observed subset with closely matching Mb and iden- 
tical T-type (we use three Mb bins for each T-type; the 
selected galaxy can be any of the observed galaxies in same 
bin) and using the corresponding SFR. 

This method requires the parent and observed samples 
to have well-matched properties at each T-type, so that ev- 
ery parent galaxy can have SF information supplied from an 
observed galaxy with similar characteristics. Figure [4] shows 
the distribution of B-band absolute magnitudes for each T- 
type for both parent and observed samples, showing that 
indeed the distributions are well-matched, as would be ex- 
pected given that 40% of the full sample were observed with 
a good spread over all T-types. Kolmogorov-Smirnov (K-S) 
tests applied to each galaxy type show that the observed 
sample is consistent with being drawn randomly from the 
parent sample. 
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Faint 


Intermediate 


Bright 


T 


= 


- 2 


-9.78 


-9.63 


-10.00 


T 


= 3 


- 7 


-9.64 


-9.67 


-9.63 


T 


= 8- 


- 10 


-9.88 


-9.74 


-9.91 



Table 1. log(SFRTOT/-^B,TOT) values for observed sample galaxies, binned by luminosity and Hubble T-type. 
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Fig. 5. The Ha luminosity function derived from our 
parent sample, with type-dependent extinction correc- 
tions. The dashed line shows Ha luminosity function of 
iGallego et alj (|l995f h The horizontal scale is given both in 
terms of Log(Ha )luminosity, in Watts (bottom scale), and 
Log(SFR), in solar masses per year (top scale). 
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T— type 

Fig. 4. The distribution of absolute magnitudes amongst 
each galaxy T-type for both the parent (a) and the observed 
(b) samples. 



The following results for the total SFR density were 
obtained from the Monte Carlo analysis. Using the 
B-band luminosity function derived from the parent 
sample and type-dependent extinction corrections from 
Paper II in converting Ha fluxes to SFRs, the SFR 
density is 0.0229±0.0016 M Q yr" 1 Mpc -3 . Using the 
SDSS blue-sequence B-band luminosity function, with the 
same Ha extinction corrections, this value is reduced to 
0.0130±0.0009 Mq yr" 1 Mpc" 3 . In both cases, the scatter 
in values found using repeated realisations of the Monte 
Carlo method are less than the formal errors, which arc 
dominated by Poisson statistics on the size of the observed 
samples. It is encouraging to note that in both cases, the 
Monte Carlo method gives SFR density values within a few 
percent of those found using the SFR:Ls correlation. 



The Ha luminosity function derived from our B-band 
luminosity function, and with the type-dependent extinc- 
tion corrections, is sho wn in Fig. [5) For comp arison, the Ha 
luminosity function of IGallego et al.l (|1995f ) is overplotted 
as a dashed line, showing good agreement in overall level 
over the rang e of SFR probed bv IGallego et~all (|l995h . 

We next used the Monte Carlo technique to evalu- 
ate the effect of different Ha extinction corrections on 
the derived SFR density, using the -B-band luminosity 
function from the parent sample of the present study. 
Replacing the type-dependent co rrection wi t h the constant 
1.1 mag correction proposed in iKennicuttl (|l998l ) gives a 
SFR density of 0.0203±0.0014 M Q yr" 1 Mpc" 3 . A third 
op tion is the Mjf-depende nt extinction correction derived 
bv lHelmboldt et all (|2004f ) based on the study of a 21 cm- 
selected galaxy sample (we also use their [Nil] contamina- 
tion correction, for consistency). This gives rather lower 
SFRs for a given Ha flux, and a total SFR density of 
0.0159±0.0012 M Q yr" 1 Mpc" 3 . Thus the systematic un- 
certainties associated with internal extinction corrections 
are ±20%, and much larger than the formal statistical er- 
rors. 
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Fig. 6. SFR density as a function of absolute blue magni- 
tude found using the points in the £?-band luminosity func- 
tion shown in Fig. [2] (top); and using the SDSS £?-band 
luminosity function (bottom). The cumulative fraction of 
the SFR density is shown by the dot-dashed line in each 
plot; the scale for this fraction is marked along the right- 
hand axis. 



5. The contribution of galaxies of different types 
and luminosities to the local SFR density 

5.1. SFR density as a function of galaxy luminosity 

Using the Monte Carlo method described in the previous 
Section, it is simple to evaluate the contributions to the 
local SFR density from galaxies of different absolute mag- 
nitudes. This is of interest given the relatively large number 
of dwarf galaxies in the present sample, which could poten- 
tially offset the low individual SFRs of these dwarfs. 

This calculation was done using both the f?-band lumi- 
nosity function derived from our parent sample, and that 
of the SDSS blue sequence. The results are shown in Fig. [6l 
from which it is clear that the local SFR density is domi- 
nated by bright spiral galaxies with Mb ~ —19 — 20 mag, 
regardless of which luminosity function is used. Even with 
our relatively dwarf-rich luminosity function, only 8% of 
the SF in the local Universe is occurring in dwarf galaxies 
fainter than Mb =-15.5 mag; using the SDSS blue-sequence 
luminosity function this falls to 6%. 




T type 

Fig. 7. SFR density as a function of galaxy morphological 
type, with galaxy numbers drawn from our i?-band lumi- 
nosity function. 

5.2. SFR density as a function of galaxy type 

We can also use our Monte Carlo method to explore the 
fraction of the total SFR density contributed by galaxies 
of different types. This uses a similar method to that in 
the previous Section (5.1), but makes use of the T-type 
information that is available for all galaxies in the parent 
sample used in this simulation. 

The results are shown in Fig. which shows that the 
galaxies contributing the most to the local SF are those 
classified as Sb - Sc. Early type spirals provide only a small 
contribution in the present epoch, in spite of including 
many of the highest mass spiral galaxies. It is interesting 
to note the large step in SFR density between T —2 and 3 
(Sab and Sb), since these types also straddle the transition 
between the SDSS red and blue sequences, as explained in 
Appendix A. Late-type Magellanic galaxies of types Sm & 
Im (T = 9 and 10), which are the most numerous galax- 
ies in the field, contribute ~13% (our £?-band luminosity 
function) or ~10% (SDSS luminosity function) of the local 
SFR density. 



6. Disk vs central SF 

A final study, made possible by the spatial coverage and 
resolution of the Ha imaging technique, is to separate out 
SF taking place in central and disk regions of the galax- 
ies observed. In this first analysis, we do not attempt a 
bulge-disk deconvolution of the galaxies, but simply define 
the central region ('bu lge') as being the central 1 kpc in 
radius (following, e.g., iGalaz et al.ll2006f ). with 'disk' SF 
being traced by the Ha emission lying outside this radius. 
For inclined galaxies, elliptical apertures with semi-major 
axis 1 kpc and axial ratios matching those of the outer disk 
were used, as discussed in Paper I. Disk and bulge fractions 
of SF were calculated in this way for all galaxies of types 
SOa - Sdm in our observed sample; the latest-type galaxies 
(Sdm, Sm and Im) are deemed to be bulge-free, and hence 
the disk fraction for these galaxies is 100% by definition. 

One reason for doing this analysis is that the disk emis- 
sion thus defined can confidently be interpreted as being 
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dominated by SF. The bulge component can potentially be 
contaminated by AGN-related emission, and this process 
puts an upper limit on the potential over-estimate of our 
SFR density values from this cause. Similarly, central re- 
gions can be affected by high [Nil] /Ha ratios (Paper II and 
references therein), which can again lead to over-estimates 
of SFRs. On the other hand, by interpreting all the emis- 
sion (bulge and disk) as being due to Ha from SF, we can 
derive an estimate of the relative rates of growth of the 
stellar mass of the two components in the current epoch, 
bearing in mind that the bulge estimates are likely to be 
upper limits. 

We find the following disk percentages of the local SFR 
density as a function of galaxy T-type: T=0, 63%; T=l, 
54%; T=2, 73%; T=3, 64%; T=4, 76%; T=5, 83%; T=6, 
81%; T=7, 78%. Types T=8, 9 and 10 are defined to have 
disk SF fractions of 100% (pure disk systems). The total 
fraction of the SFR density occurring within disks is 75% 
for types 0-7; this increases to 80% when including the SF 
of the later type galaxies in the disk total. This demon- 
strates that the major part of the SF activity is occurring 
in disk regions, and hence that our estimates of the total 
SFR density cannot be significantly affected by an AGN 
contribution. The bulge fraction of line emission for the 
earliest-type galaxies is high, however, at almost 50% of the 
total. This may imply that these galaxies are still building 
the stellar mass in their central regions at a significant rate, 
or it may be the result of substantial contamination from 
AGN or other non-SF related emission. 

A fuller investigation of the spatial distribution of line 
emission within these galaxies will be described in a later 
paper in this series (James et al., in preparation). 

7. Comparison with other determinations of SFR 
density 

We finally compare our estimates of the local SFR den- 
sity with estimates of this same parameter from the lit- 
erature, where the latter studies use a variety of field 
galaxy samples and tracers of galaxy SFRs. The present 
study found values of this parameter ranging from 0.0159 
to 0.0229 Mq yr -1 Mpc -3 using the internally-derived B- 
band luminosity function, with the range reflecting the dif- 
ferent internal extinction corrections used. Using the B- 
band luminosity function for the SDSS blue sequence gave 
a lower estimate, of 0.0130 M Q yr _1 Mpc -3 . 

Liter ature results in this area have recently been sur- 
veyed by lHanish et al.1 (|2006l ). who quote results from sev- 
eral surveys, with and without extinction corrections, and 
with and without a redshift correction, since the surveys 
quoted have mean effective redshifts from 0.01 to 0.20. The 
numbers presented here include both extinction correction 
and the redshift correction to bring the SFR densities to an 
equivalent at z = 0. 

Two results come from the Universidad Complutense 
de Madrid su rvey of opti c ally-se lected emission line 
galaxies, with iGallego et alj (|1995l) finding a value of 
0112 ±0.001 M Q yr" 1 Mpc" 3 , and iPerez-Gonzalez et all 
(12003ft 0.0229±0.001 M Q yr" 1 M pc" 3 . The /-band s e lected 
Canada- France Redshift survey of lTresse fc Maddoxl fl998) 
resulted in a SFR density value of 0.0186±0.0032 M Q 
yr -1 Mpc -3 , co rrected from an effec tive mean redshift 
of 0.2 to z = 0. ISullivan et all (|200Ct ) use a UV-selected 



sample of galaxies over quite a wide redshift range (0 
- 0.4) to derive a SFR density of 0.01±0.001 M© yr" 1 
Mpc -3 , although this low value is comp letely depen- 
dent on the redshift correction adopted by lHanish et al.l 
(2006). Easily the most extensive study in terms of num- 
bers of galaxies in clude d, is the SDSS-based study of 
iBrinchmann et all (|2004l ). who find a SFR density of 
0.0219^^048 M Q vr " 1 Mpc" 3 . Finally for their own, Hl- 
selected field sample, lHanish et al.1 (|2006l ) find a SFR den- 
sity of 0.0155±0.001 M Q yr -1 Mpc -3 . Overall, the range 
of values found is in good agreement with those derived in 
the current paper, and there is no clear reason to favour 
either the high values we obtain with the HaGS luminosity 
function, or the lower SDSS-derived value. 

8. Summary 

We have presented an analysis of the SF properties of field 
galaxies within the local volume out to a recession velocity 
limit of 3000 km s _1 . A parent sample of 863 galaxies with 
well-defined Hubble type, apparent magnitude and diame- 
ter limits was selected from the UGC, and used to calcu- 
late a i3-band luminosity function. This was then populated 
with SF information from a subsample of 327 galaxies, for 
which we have Ha imaging, firstly by calibrating a relation- 
ship between galaxy _B-band luminosity and SFR, and sec- 
ondly by a Monte Carlo simulation of a representative sam- 
ple of galaxies, in which SF information was randomly sam- 
pled from the observed subset. The total SFR density of the 
local Universe was found to be between 0.016 and 0.023 Mq 
yr -1 Mpc~ 3 regardless of which of these two methods was 
used, with the uncertainties being dominated by the in- 
ternal extinction correction used in converting measured 
Ha fluxes to SF rates. If our internally derived B-band lu- 
minosity function is replaced by one from the SDSS blue 
sequence, the SFR densities are systematically reduced to 
about 60% of the above values. This range of estimates is 
consistent with values from the recent literature using a 
range of different SF indicators. 

The Monte Carlo method was then used to calculate the 
contribution to the total SFR density from galaxies of dif- 
ferent luminosities and Hubble T-types. The largest contri- 
bution comes from bright galaxies with Mb ^-20 mag, and 
the total contribution from galaxies fainter than Mb = 
15.5 mag is less than 10%. Almost 60% of the SFR den- 
sity comes from galaxies of types Sb, Sbc or Sc; 9% from 
galaxies earlier than Sb and 33% from galaxies later than 
Sc. Finally, 75 - 80% of the total SF in the local Universe 
is shown to be occurring in disk regions, defined as being 
>1 kpc from the centres of galaxies. 
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Appendix A: Morphological classification along the 
red and blue sequences 

With the advent of large imaging surveys, both local and 
at high redshift, colour classification in particular dividing 
into a red and blue sequence is often preferred to morpho- 
logical classification. The is partly because of the reliability 
of determining galaxy colours automatically compared to 
morphological characteristics for th ese surveys, and partl y 
because of the colour bimodality (Str ateva et all [2001). 
Galaxies divide naturally into a red and blue sequence when 
viewed in a colour-mag nitude diagram (CMP) using u — r , 
q — r colours or similar (jBlanton et all 12003: Bal dry et all 
2004). In order to outline the connection between tradi- 



Table A.l. Fraction of blue-sequence galaxies for different 
morphological types 



morphological 


NFGS 




SDSS JPG 


group 


no. 


h 


no. 


A (A) 


A(B) 


E 


17 


0.24 


201 


0.04 


0.03 


SO 


28 


0.21 


114 


0.18 


0.11 


SO/a 


12 


0.42 


73 


0.19 


0.22 


Sa,ab 


21 


0.52 


131 


0.35 


0.37 


Sb,bc 


31 


0.81 


152 


0.61 


0.73 


Sc,cd 


45 


0.93 


121 


0.86 


0.92 


Sd,dm 


18 


0.94 


12 


0.83 


0.75 


Im 


22 


1.00 


10 


0.90 


1.00 


all 


194 


0.68 


814 


0.37 


0.40 



(A) using cut in [u - r) mode i CMD, (B) using cut in (g - r) PctI , 
CMD. 



tional morphological classification using the E-SO-Sa-Sd-Im 
scheme and the colour sequences, we used data from the 
Nearby Field Galaxy Survey (NFGS; I Jansen et all l2000h 
and t he SDSS Japanese Pa rticipation Group (JPG) cata- 
logue d^kugitaiiin,[l003). 

The NFGS includes UBR surface photometry for 198 
galaxies with a median redshift of 0. 01. The NFGS was 
selected from the iHuchra et al.l f|l983f ) catalogue, exclud- 
ing the Virgo cluster, with sampling to obtain a maxi- 
mum of about 30 galaxies per 1-mag bin (blue absolute 
magnitude). We used the U — R colours within the effec- 
tive half-light radius in B, converted to AB magnitudes 
(+0.55), and used the absolute £?-band magnitude con- 
verted to H = 70 km s _1 Mpc" 1 . Fi gure IA.f1 (left) shows 
the CMD for these data with different symbols representing 
morphological classification in eight groups. 

The JPG catalogue includes SDSS photometry and eye- 
ball classifications for 2658 galaxies. This was reduced to 
818 galaxies with r < 15.5 mag and z < 0.1 (median red- 
shift of 0.04). We used the u — r 'model' rest-frame colour 
and a .B-band absolute magnitude estimated using the g 
and r magnitudes. SDSS mo del magnitudes are bas ed on a 
fitted profile in the r-band (|Stoughton et all l2002f) : effec- 
tively they are centrally weighted colours. Figure [A~lT right) 
shows the CMD for these data after sampling to obtain a 
maximum of 25 galaxies per 0.5-mag bin. The data are di- 
vided into eight classification groups (note that an interme- 
diate classification, e.g. Sab, represents an average between 
two or more classifiers.) 

The dashed line in Figure IA.1I re presen ts the best-fit 
dividing line defined by iBaldrv etal] |2006l ). This was de- 
termined from a significantly larger SDSS sample. This ap- 
pears to fall in about the right place for the NFGS despite 
the different colour determination. The faint-end of the red 
sequence is sparsely populated because these are field sam- 
ples. In order to compare late- type and blue-sequence de- 
rived luminosity functions, we determined the fraction of 
galaxies on the blue sequence (/&) for the eight morpholog- 
ical groups. This is shown in Table IA~T1 

When comparing fa it should first be noted that the 
NFGS and JPG samples were selected using different cri- 
teria. For example, the JPG sample is r-band magnitude 
limited and there was no explicit removal of galaxies in 
high-density regions. Therefore, it is not surprising that the 
overall is lower. Another factor is the higher median red- 
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Fig. A.l. Morphological classifications in the CMD for galaxies. Left: NFGS sample. Right: A subset of the JPG sample. 
The dashed line in both plots represents an optimal divider between the red and blue sequences determined from a larger 
SDSS sample. 



shift of the JPG which affects the quality of morphological 
classification (at moderate imaging resolution of ~ 2"). For 
the purposes of this paper the most important difference is 
the fb fraction for the Sb,bc galaxies because of their high 
contribution to the SFR density, as found in section 5.2. 
This value is 0.61 for the JPG sample (compared to 0.81 
for the NFGS). This may partly account for the lower SFR 
density derived from the blue-sequence SDSS luminosity 
function (section 4). 

One possible reason for a low fb value for the Sb- 
classified JPG galaxies is the use of the u — r model colour, 
which is centrally weighted. It is generally not reliable to 
use SDSS u-band magnitudes derived from larger apertures 
because of Poisson noise and sky subtraction, and an alter- 
native is the g — r Petrosian colour. Figure [A~2l shows the 
g — r CMD for the JPG sample, and the fb values are also 
shown in the final column of Table The value is now 
0.73 for Sb,bc galaxies indicating that aperture choice not 
surprisingly affects the red and blue sequence classification 
of these galaxies. Nevertheless, a blue-sequence luminosity 
function can be regarded as being similar to an Sb-or-later 
luminosity function because of the cross-over in fb between 
Sa and Sb (considering the magnitude-limited JPG sam- 
pie). 



12 James et al.: Ha Galaxy Survey. IV. 



1.0 = 



+ + O A V □ X 
E SO SO/a Sa Sb Sc Sd 
E.. Sab Sbc Scd Sdm 



X 



- -B - 



- - - e + ^ v ™ qa <y □ + 



< 



Fig. A. 2. Morphological classifications in the CMD for galaxies. As per Figure [A~T1 (right) except using g — r Petrosian 
colours instead of u — r model colours. The dashed line represents an optimal divider determined by eye from a larger 
SDSS sample: (g - r) dividc = 0.65 - 0.023(M b + 20). 



